Decoherence in Attosecond Photoionization 
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The creation of superpositions of hole states via single-photon ionization using attosecond extreme- 
ultraviolet pulses is studied with the time- dependent configuration interaction singles (TDCIS) 
method. Specifically, the degree of coherence between hole states in atomic xenon is investigated. 
We find that interchannel coupling not only affects the hole populations, it also enhances the entan- 
glement between the photoelectron and the remaining ion, thereby reducing the coherence within 
the ion. As a consequence, even if the spectral bandwidth of the ionizing pulse exceeds the energy 
splittings among the hole states involved, perfectly coherent hole wave packets cannot be formed. 
For sufficiently large spectral bandwidth, the coherence can only be increased by increasing the 
mean photon energy. 

PACS numbers: 32.80.Aa, 42.65.Re, 03.65.Yz 



The typical time scale of electronic motion in atoms, 
molecules, and condensed matter systems ranges from a 
few attoseconds (1 as = 10""'^* s) to tens of femtoseconds 

(1 fs = 10-15 s) dm. In the last decade the remark- 
able progress in high harmonic generation made it 
possible to generate attosecond pulses as short as 80 as 
[9| . Attosecond pulses have opened the door to real-time 
observations of the most fundamental processes on the 
atomic scale [3, [l^l . For instance, the generation of at- 
tosecond pulses was utilized to determine spatial struc- 
tures of molecular orbitals [11]; an intcrfcrometric tech- 
nique using attosecond pulses was used to characterize 
attosecond electron wave packets [1^; and attosecond 
pulse trains [l3| and isolated attosecond pulses il4], in 
combination with an intense few-cycle infrared pulse, en- 
abled the control of electron localization in molecules. 
Attosecond technology demonstrated the ability to fol- 
low, on a subfemtosecond time scale, processes such as 
photoionization jT^, Auger decay [l6|, and valence elec- 
tron motion driven by relativistic spin-orbit coupling [l7| . 
Furthermore, the availability of attosecond pulses fuelled 
a broad interest in exploring charge transfe r dy namics 
following photoexcitation or photoionization 14| . 



In this Letter, we analyze the creation of hole states 
via single-photon ionization using a single extreme- 
ultraviolet attosecond pulse. We investigate the impact 
of the freed photoelectron on the remaining ion and 
demonstrate that the interaction between the photoelec- 
tron and the ion cannot be neglected for currently avail- 
able state-of-the-art attosecond pulses. In particular, the 
interchannel coupling of the initially coherently excited 
hole states greatly enhances the entanglement between 
the photoelectron and the ionic states. Interchannel cou- 
pling is mediated by the photoelectron and mixes dif- 
ferent ionization channels, i.e., hole configurations, with 



each other. Consequently, the degree of coherence among 
the ionic states is strongly reduced, making it impossible 
to describe the subsequent charge transfer in the ion with 
a pure quantum mechanical state. Experiments on pho- 
tosynthetic systems 18l-l2l| have revealed a correlation 
between highly efficient energy transport and coherent 
dynamics in molecules (nuclear and electronic dynam- 
ics in this case). Similarly, high degrees of coherence in 
nonstationary hole states may be necessary for efficient 
charge transport within molecules. 

In the last decade, much work has been done in the 
realm of hole migration [22 -24]. It was shown that elec- 
tronic motion can be triggered solely by electron corre- 
lation [22I . Charge transfers mediated by electronic cor- 
relations typically take place in a few femtoseconds and 
are thus faster than electronic dynamics initiated by nu- 
clear motion [111, [26| . Recent experiments 27 , [i^] have 
demonstrated that electronically excited ionic states can 
modify site-selective reactivity within tens of femtosec- 
onds, making hole migration processes a promising tool 
to control chemical reactions. Up to now, theoretical cal- 
culations 22 , 23] investigating hole migration phenomena 
have neglected the interaction between the parent ion and 
the photoelectron and assumed a perfectly coherent hole 
wave packet. As long as the photoelectron departs suf- 
ficiently rapidly from the parent ion, this assumption is 
appropriate [29| . However, for attosecond pulses with 
large spectral bandwidths, the enhanced production of 
slow photoelectrons will affect (mainly via interchannel 
coupling) both the final hole populations and the coher- 
ence among these hole states. Furthermore, recent results 
in high harmonic spectroscopy suggest that interchannel 
coupling may be the missing link to understand hole dy- 
namics occurring in high harmonic generation processes 
before the ejected electron recombines with the parent 
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ion [30|. 

We investigate the creation of hole states via attosec- 
ond photoionization using the implementation of the 
time-dependent configuration-interaction singles (TD- 
CIS) approach described in Ref. 



31| (see also 
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TDCIS allows us to study ionization dynamics beyond 
the single-channel approximation and to understand sys- 
tematically the relevance of interchannel coupling in the 
hole creation process. The TDCIS wave function for the 
entire system is 



\^it)) = aoit) |<fo) 



(1) 



where |$o) is the Hartree-Fock ground state and |$°) = 
c^Ci |$o) is a one-particle-one- hole excitation (cj^ and Ci 
are creation and annihilation operators for an electron 
in orbitals a and i, respectively). The corresponding co- 
efficients ao{t) and af{t)^ respectively, are functions of 
time and describe the dynamics of the system. Through- 
out, indices are used for occupied orbitals in |$o); 
indices a, 6, stand for unoccupied orbitals. We focus our 
discussion on the case where single-photon ionization is 
the dominant effect and higher order processes can be 
neglected. Our model system is atomic xenon. The cor- 
responding Hamiltonian (neglecting spin-orbit coupling) 
is 



H{t) = Ho + Hi+E{t)^ 



(2a) 



where E{t) is the electric field, z the dipole operator, and 
Ha is the mean-field Fock operator, which is diagonal 
with respect to the basis used in Eq. ([T]). The residual 
Coulomb interaction. 



Hi — Vc — Vmf, 



(2b) 



is defined such that Hq + Hi gives the exact nonrelativis- 
tic Hamiltonian for the electronic system in the absence 
of external fields {Vc is the electron-electron interaction). 
We study the impact of different approximations for Hi 
on the hole state as follows. The Coulomb-free model, the 
simplest approximation, removes the residual Coulomb 
interaction {Hi = 0) between the excited electron and 
the parent ion. In this approximation, the excited elec- 
tron always sees a neutral atom via the Vmf potential 
34| . A more realistic approximation is the intrachan- 
nel model including direct and exchange contributions of 
the Coulomb interaction only within a given channel. In 
this second model, the excited electron can only inter- 
act with the occupied orbital from which it originates. 
Interactions between different occupied orbitals are ne- 
glected, i.e. we set j^J) = for i ^ j. The 
third and final model describes the Coulomb interaction 
exactly within the TDCIS framework. We refer to this 
as the full model. Note that the exact nonrelativistic 
Hamiltonian Hq + Hi is diagonal with respect to the 



ionic one-hole states = Ci |$o)- In the full model, 
the photoelectron can couple the hole states, as Hi in 
the particle-hole space is not diagonal with respect to 
the hole index (i.e., {^f\Hi \ generally differs from 
zero). This type of photoelectron-mediated interaction 
is called interchannel coupling [sHj. As a consequence, 
in the full model the hole index is not a good quantum 
number, whereas in the Coulomb-free and intrachannel 
models, excited eigenstates of Hq + Hi are characterized 
by a well-defined hole index. To describe the hole states 
of the remaining ion, we employ the ion density matrix 



pi;j^(i) = Tr„[|*(t))(vl/(t)| 



where Tr^ stands for the trace over the photoelectron. 
The properties of the ion density matrix can be measured 
using attosecond transient absorption spectroscopy [17[. 



A description of the cationic eigenstates in terms of one- 
hole configurations is a physically meaningful approxima- 
tion for noble-gas atoms such as xenon 
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FIG. 1: (color online) The 4do [panel (a)] and 5s [panel (b)] 
hole populations of xenon as a function of time are shown 
for three different residual Coulomb interaction approxima- 
tions: (1) the full model (red solid line), (2) the intrachannel 
model (green dotted line), and (3) the Coulomb- free model 
(blue dash-dotted line). The attosecond pulse has a peak 
field strength of 25 GV/m, a pulse duration of 20 as, a (mean) 
photon energy of 136 eV, and is centered at t = as. 



In Fig. [I] the hole populations p'^,%{t) and pS„^4do(^) 
of the xenon 5s and Ado orbitals, respectively, are shown 
for all three interaction models {Ado stands for the 4c? 
orbital with m = 0). The ionizing, gaussian-shaped at- 
tosecond pulse is linearly polarized and has a peak field 
strength of 25 GV/m, a pulse duration of r = 20 as, 
and a (mean) photon energy of ujo = 136 eV. The hole 
dynamics of the Coulomb-free and intrachannel models 
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are alike. In both cases, the population is constant af- 
ter the pulse, since the hole index is a good quantum 
number within these models. The extension to the exact 
Coulomb interaction changes the situation. Interchannel 
coupling causes the hole populations to remain nonsta- 
tionary as long as the photoelectron remains close to the 
ion. As the distance between the photoelectron and the 
ion increases, the interchannel coupling weakens and the 
populations pYl^ {t) become stationary (see Fig. [TJ . We 
confine our discussion to the first hundreds of attoseconds 
after the pulse, allowing us to neglect decay processes, 
which start to take place after a few femtoseconds. 

As we will see in the following, interchannel coupling 
not only affects the hole populations but also the co- 
herence between the created hole states. The degree of 
coherence between and is given by 
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Totally incoherent statistical mixtures result in gij{t) — 
0. The fact that the density matrix is positive scmidef- 
inite implies the Cauchy-Schwarz relations ^ 
pYi^ (t) pYj^ (t) J which bound the maximum achievable 
(perfect) coherence {gi,j{t) = 1). To investigate the ef- 
fect of interchannel coupling on the coherence between 
the orbitals 4do and 5s in xenon, we restrict the defini- 
tion of the Ado hole population to the events where the 
photoelectron has angular momentum 1 = 1. The other 
possible angular momentum for the Ado photoelectron, 
Z = 3, does not contribute to the coherence, since the 
photoelectron from 5s can only have 1 = 1. For a similar 
reason, it is impossible to create a coherent superposition 
of bp and 5s (or Ad) hole states via one-photon absorption 
in the electric dipole approximation. 

Figure [2] illustrates the time evolution of the coherence 
between Ado and 5s in xenon for different pulse dura- 
tions and fixed photon energy (wq = 136 eV). Here, we 
use the full interaction model. Directly after the ion- 
izing pulse is over, the initial degree of coherence (at 
i ss as) rises with decreasing pulse duration, i.e., in- 
creasing spectral bandwidth, and converges to a value 
close to unity. (The difference of the ionization poten- 
tials, £5s — £4do, is « 50 eV.) At t « as, the pho- 
toelectron is still in immediate contact with the parent 
ion. Therefore, the coherence properties of the system 
of interest — the parent ion — are affected by its interac- 
tion with the bath represented by the photoelectron. The 
system-bath interaction leads to a reduction in the co- 
herence of the system 37[ , which can be seen by the rapid 
drops in all curves in Fig. [2] within tens of attoseconds 
after the pulse. With time, as the photoelectron departs 
from the ion, the Coulomb ("system-bath") interaction 
becomes less important and the coherence converges to a 
stationary value. The maximum for this stationary value 
is obtained with a 25 as pulse {gMo.Ss ~ 0.6). For pulses 
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FIG. 2: (color online) The time evolution of the coherence 
between the Ado and 5s hole states in xenon is shown for the 
full Coulomb interaction model. The photon energy is 136 eV 
and the pulse duration varies from 5-60 as. 



shorter than 25 as, oscillations in .94^0, 5s occur that per- 
sist for hundreds of attoseconds, and the final degree of 
coherence reached falls below 0.6. The spin-orbit dynam- 
ics associated with the fine-structure within the 4c? shell 
is slow in comparison to the time scale of the decoher- 
ence between Ado and 5s, and is, therefore, not considered 
here. 




400 600 
time [as] 

FIG. 3: (color online) The time evolution of the coherence 
between the Ado and 5s hole states, calculated with the full 
Coulomb interaction model, is shown for different photon en- 
ergies. The pulse duration is in all cases 20 as. 



We see in Fig. [3] that when holding the pulse duration 
fixed (r = 20 as), the degree of coherence rises with in- 
creasing wq. The magnitude of the oscillations decreases 
as the final coherence (at i w 1 fs) increases. This trend 
indicates less system-bath interactions occur with higher 
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photoelectron energies keeping the degree of coherence 
among the hole states high. 

In Fig. m we compare the impact of the different 
Coulomb approximations on the final coherence. The 
drops in coherence that occur for the full model for short 
pulses [Fig. mja)] and low photon energies [Fig. Hl^b)] 
cannot be seen in the Coulomb-free and intrachan- 
nel models — which both neglect interchannel coupling. 
Hence, the decay of coherence is solely driven by the in- 
terchannel coupling due to the slow photoelectron. As 
a comparison to the Coulomb-free model shows, intra- 
channel coupling affects the coherence in an insignificant 
way. In the limit of long pulse durations (small spectral 
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FIG. 4: (color online) The dependence of the coherence be- 
tween the 4do and 5s hole states as function of the pulse du- 
ration (a) and as function of the photon energy (b) are shown 
for all three interaction approximations. 

bandwidths), the coherence vanishes for all models, since 
photoelectrons from the Mq and 5s become energetically 
distinguishable and cannot contribute to a coherent sta- 
tistical mixture of hole states. The slight drop in the 
coherence for the Coulomb-free and intrachannel models 
with increasing wq [Fig. IHJb)] is related to the reduced 
factorizability of the numerator of Eq. In contrast, 
the trend in the full model for increasing luq is dominated 
by the gain in coherence due to higher photoelectron en- 
ergy resulting in less system-bath interaction. 

In conclusion, we demonstrated that the coherence of 
the ionic states produced via attosecond photoionization 
is not solely determined by the bandwidth of the ionizing 
pulse, but greatly depends on the kinetic energy of the 
photoelectron, which can be controlled by the (mean) 
photon energy. Interchannel coupling leads to an en- 
hanced entanglement between the photoelectron and the 
parent ion resulting in a reduced coherence in the ionic 
states. This reduction can be mitigated with higher pho- 
ton energies, thereby sacrificing high photon cross sec- 
tions and the possibility of controlling independently the 
relative populations of the various hole states in the sta- 
tistical mixture. 



Our results have far-reaching consequences beyond the 
atomic case. Molecules will be even more strongly af- 
fected by interchannel coupling due to the reduced sym- 
metry and smaller energy splittings between the cation 
many-electron eigenstates. Interchannel coupling is also 
likely to be significant for inner-valence hole configura- 
tions in molecules, which show strong mixing to config- 
urations outside the TDCIS model space. The present 
study suggests that interchannel coupling accompany- 
ing the hole creation process will affect attosecond ex- 
periments investigating charge transfer processes in pho- 
toionized systems. The control of decoherence requires 
widely tunable attosecond sources, thus offering a new 
opportunity for x-ray free-electron lasers [s^j. 

P.J.H. was supported by the Office of Basic Energy 
Sciences, U.S. Department of Energy under Contract 
No. DE-AC02-06CH11357. L.G. thanks Martha Ann 
and Joseph A. Chenicek and their family. D.A.M. grate- 
fully acknowledges the NSF, the Henry-Camille Dreyfus 
Foundation, the David-Lucile Packard Foundation, and 
the Microsoft Corporation for their support. 



* Corresponding author; Electronic address: 

robin.santra@cfel.de 
[1] F. Krausz and M. Ivanov, Rev. Mod. Phys. 81, 163 
(2009). 

[2] A. H. Zewail, J. Phys. Chem. A 104, 5660 (2000). 
[3] A. L. Cavalieri et ai, Nature 449, 1029 (2007). 
[4] G. Doumy et al, Phys. Rev. Lett. 102, 093002 (2009). 
[5] N. Dudovich et al, Nature Phys. 2, 781 (2006). 
[6] R. Lopez-Martens et ai, Phys. Rev. Lett. 94, 033001 
(2005). 

[7] E. A. Gibson et al, Phys. Rev. Lett. 92, 033001 (2004). 

[8] K. J. Schafer et ai, Phys. Rev. Lett. 70, 1599 (1993). 

[9] E. Goulielmakis et al, Science 320, 1614 (2008). 
[10] T. Pfeifer et al, Chem. Phys. Lett. 463, 11 (2008). 
[11] S. Haessler et al, Nature Phys 6, 200 (2010). 
[12] J. Mauritsson et al, Phys. Rev. Lett. 105, 053001 (2010). 
[13] K. P. Singh et al, Phys. Rev. Lett. 104, 023001 (2010). 
[14] G. Sansone et al., Nature 465, 763 (2010). 
[15] E. Skantzakis et al, Phys. Rev. Lett. 105, 043902 (2010). 
[16] M. Drescher et al. Nature 419, 803 (2002). 
[17] E. Goulielmakis et ai. Nature 466, 739 (2010). 
[18] H. Lee, Y.-C. Cheng, and G. R. Fleming, Science 316, 
1462 (2007). 

[19] M. Sarovar et al. Nature Phys. 6, 462 (2010). 

[20] E. CoUini et ai. Nature 463,' 644 (2010). 

[21] E. Harel, A. F. Fidler, and G. S. Engel, Proc. Natl. Acad. 

Sci. USA 107, 16444 (2010). 
[22] J. Breidbach and L. S. Cederbaum, J. Chem. Phys. 118, 

3983 (2003). 

[23] M. Nest, F. Remade, and R. D. Levine, New J. Phys. 

10, 025019 (2008). 
[24] A. Kuleff, S. Liinnemann, and L. Cederbaum, J. Phys. 

Chem. A 114, 8676 (2010). 
[25] D. Kroner et al, Appl. Phys. A-Mater. 88, 535 (2007). 
[26] B. H. Muskatel, F. Remade, and R. D. Levine, Phys. 



5 



Scripta 80, 048101 (2009). 
[27] R. Wcinkauf et al., J. Phys. Chem. A 101, 7702 (1997). 
[28] E. W. Schlag et al., Angew. Chem. Int. Ed. 46, 3196 

(2007). 

[29] L. Cederbaum et al., Adv. Chem. Phys. 65, 115 (1986). 
[30] Y. Mairesse et al, Phys. Rev. Lett. 104, 213601 (2010). 
[31] L. Greenman et al., Phys. Rev. A 82, 023406 (2010). 
[32] H. B. Schlegel, S. M. Smith, and X. Li, J. Chem. Phys. 

126, 244110 (2007). 
[33] S. KUnkusch, P. Saalfrank, and T. Klamroth, J. Chem. 

Phys. 131, 114304 (2009). 



[34] A. Szabo and N. S. Ostlund, Modem Quantum Chemistry 

(Dover, Mincola, NY, 1996). 
[35] A. F. Staracc, Appl. Opt. 19, 4051 (1980). 
[36] C. Buth, R. Santra, and L. S. Cederbaum, J. Chem. Phys. 

119, 7763 (2003). 
[37] H.-P. Breuer and F. Petruccione, The Theory of Open 

Quantum Systems (Oxford University Press, 2002). 
[38] A. Zholents and W. Fawley, Phys. Rev. Lett. 92, 224801 

(2004). 



